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Abstract 

Based on a proposed classical explanation, the quantum mechanical "decay of the wave packet" is 
shown to simply result from sub-quantum diffusion with a specific diffusivity varying in time due to 
a particle's changing thermal environment. The exact quantum mechanical intensity distribution, 
as well as the corresponding trajectory distribution and the velocity field of a Gaussian wave packet 
are therewith computed. We utilize no quantum mechanics, but only familiar simulation techniques 
for diffusion, e.g., finite differences or coupled map lattices (CML). 
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1. BALLISTIC DIFFUSION 



We derive a solution for an anomalous diffusion equation with a time-dependent diffusion 
coefficient D t (t) = kt a of the general form 

OP d 2 P 

Here, t and /c denote the time and a constant factor, respectively. P(x,t) is the solution of 
Eq. (1.1). We assume a Gaussian distribution function as ansatz for P, 

(x - x ) 2 

P(x,f) = -_e 2a 2 , (1.2) 
y 2ti<j 



where a(t) denotes the standard deviation, i.e. the distance from the mean value to the 
left or right interception point of the distribution function P. For further considerations we 
put the initial location x in Eq. (1.2) into the coordinate's origin, i.e., Xq = 0. Then we 
calculate the partial derivatives of P using the short form & for the time-derivative of cr(t), 

then substitute these results into Eq. (1.1), 

Pa fx 2 \ , „, P ( x 2 



K<J 2 j a 2 \a 2 

and obtain 

T = *STT + f- ^ 

At this stage we introduce a result of [1] which is an expression for the spreading of the 
Gaussian, 

- 2 = -S(l + 5 ! ) , (1-6) 

with o"o as initial standard deviation at time t — 0. Note that the diffusivity D = -^ L - 
is constant for all times t and has to be distinguished from the diffusion coefficient D t . 
Substitution of Eq. (1.6) into (1.5) yields c = a 2 and 

, 2t a+1 D 2 t 2 , s 

k - = —-. 1.7 

a + 1 a 2 K ' 
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Eq. (1.7) can only be fulfilled by a = 1, so that k = ^ . Hence, our time-dependent diffusion 
coefficient becomes 

A = —a- • (1.8) 



^0 



Finally, Eq. (1.1) becomes 



dP <9 2 P , s 



and turns out to be a ballistic diffusion equation because a — 1. 



2. CLASSICAL SIMULATION 

It is straightforward to simulate the diffusion process of Eq. (1.9) in a computer model. 
One approximates the diffusion equation by 

P[x, t+l] = P[x, t] + ^3 {P[x + 1, *] - 2P[g, t] + P[x - 1, t}} (2.1) 

with space and time grid indices x and t, respectively, and an initial Gaussian P(x, 0) at 
t — with initial standard deviation cr - 



3. DOUBLE SLIT INTERFERENCE 

The simulation of interference of two beams emerging from Gaussian slits is established 
with the aid of a simulation as in Refs. [1-4]. To account for interference, we simply follow 
the classical rule for the intensities P tot := R 2 = \Riki + -R2k 2 1 2 = R\ + R\ + 2RiR 2 cosip = 
Pi + P2 + 1\JP\P 2 costp, with Rk = -Rikx + R 2 ^2 an d <p = mAv x x/H . The trajectories are 
the flux lines obtained by choosing a set of initial points at y = 0. Two adjacent flux lines 
thereby define regions of constant flux, i.e., J A P dA = const., with A being the cross section 
of a flux tube. 

We have shown in [3, 4] that it is possible to simulate a particle's trajectories behind a 
double slit with usual classical simulation tools. The key is the ballistic diffusion equation 
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(1.9) which describes the underlying physics by the time-dependent diffusion coefficient (1.8). 
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